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Quantum fluctuation of unstable modes about gravitational instantons causes the insta- 
bility of flat space at finite temperature, leading to the spontaneous process of nucleating 
quantum black holes. The energy-density of quantum black holes, depending on the initial 
temperature, gives the cosmological term, which naturally accounts for the inflationary phase 
of Early Universe. The reheating phase is attributed to the Hawking radiation and annihi- 

£N| . lation of these quantum black holes. Then, the radiation energy-density dominates over the 

U 

' In this phase the energy-density of quantum black holes depends on the reheating tempera- 

ture. It asymptotically approaches to the cosmological constant in matter domination phase, 

fNl • consistently with current observations. 
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I. INTRODUCTION. 

Oh, In recent years, as observational data concerning on cosmology are rapidly accurate, the theo- 

D . n 

retical understanding of our Universe has been greatly profounded ever before. The inflation [l[ 



in early Universe and the acceleration of present Universe 2j are two most important issues in 



modern cosmology and fundamental physics. Both issues are closely related to the cosmological 
term in Einstein equation. The observational value of the cosmological constant 

A obs = 8vrGe A ~ W~ U8 m 2 p , (1) 

where the Newton constant G = 1/nip, the Planck mass m p ~ 10 19 GeV and a related characteristic 
energy-density 

6 A ~ 10- 120 mJ. (2) 

This energy-density is of 120 order smaller than the Planck energy-density, epi anc k = nip. 

Much effort in understanding these issues has been made for many decades and there are 
many interesting ideas and innovative theoretical developments based on either simple models Q] 
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or complex theories [-J. In these approaches, a scalar field slowly-rolling downwards an effective 
potential, which mimics a possible vacuum-energy variation, plays a crucial role in driving inflation; 
whereas this scalar field with a mass of order the current Hubble scale, can possibly account for 
acceleration. 

In this article, we attempt to study these two issues within a framework based on the instability 
of flat space at finite temperature, due to an unstable quantum fluctuation about gravitational 
instantons. We shall use the signature ( — h ++). 



II. COSMOLOGICAL TERM AND ENERGY-MOMENTUM TENSORS 
A. Cosmological constant term 

The gravitational field interacts only with the total energy-momentum T ab that is attributed 
to both fundamental particles and vacuum fluctuations. The total energy- momentum tensor Tab 
appears as source of the gravitational field. The gravitational field equation should therefore be of 
form 

Sab = 8TrGT ab , (3) 

where G is the Newton constant. In the left-handed side of the field equation, Gab is some expression 
involving the metric g a b and its derivatives. The field equation ([3]) is invariant under general 
coordinate transformations, therefore Gab must be a tensor. It is further assumed that the field 
equation ([3]) is of second differential order, and is linear in the second derivatives. We know that the 
only tensor that contains the second derivative linearly is R a bcd- Hence the only tensors available 
for the construction of Gab are R a bcd and g a b- The most general second-rank tensor that can be 
built out of these, with R a bcd entering linearly, is 

Gab = Rab + Bg ab R + Ag ab , (4) 

where R ab (R) is the Ricci tensor (scalar), B and A are constants. The conservation law for the 
energy-momentum tensor of matter fields T^. b = requires the object Gab must satisfy the identity 

Ga;b = 0. (5) 

In view of g a . b = and the Bianchi identity 

(R b a ~ \t% = 0, (6) 
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the identity ([5]) is fulfilled if B = 1/2. The A remains as an arbitrary constant. Thus, gravitational 
field equation becomes 

Rab ~ ^9abR + &g ab = 8irGT ab , (7) 

where Ag ab is the cosmological term (A-term) and A is the so-called cosmological constant. The 
A-term clearly has its gravitational origin (or geometric origin), rather than material origin. In 
Newtonian gravitating systems, the cosmological constant A = consistently with observations 
and experiments. 



B. Energy-momentum tensors of fundamental particles and vacuum 

In the context of local quantum field theories, the total energy-momentum tensor T ab , the 
right-handed side of the gravitational field equation (|3I7|) . can be written as, 

% b = <0| (f ab + f$) |0> + (0|Af^|0>, (8) 

where |0) stands for the vacuum state and (0| • • ■ |0) denotes the vacuum expectational value. T ab 
represents the energy-momentum of fundamental particles (including gravitational field) and T^ b 
represents the energy-momentum of the vacuum. AT^ represents the energy-momentum contri- 
bution from particles back interacting on the vacuum. 

Only energy-momentum difference between quantum states can be physically significant in the 
flat space time. As a consequence, the vacuum is defined as the ground state with zero energy- 
momentum, 

(0|T^[0) ^0 (9) 

and fundamental particles are energy-momentum excitations from the vacuum. The energy- 
momentum of fundamental particles is then given by 

T ab = (0| (f ab + fl) |0) = <0|T o6 |0>. (10) 

The T ab is in fact the net contribution from fundamental particles and the vacuum in the context of 
quantum field theories. It corresponds the classical energy-momentum, entering the right-handed 
side of the field equation (|7|8p as the source for gravitational field. 

However, we expect to have some problems in such a scenario: (i) the ground state (vacuum) of 
quantum gravity at the Planck scale is still an open question; (ii) the energy-momentum difference 
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([9]) between fundamental particles states and vacuum state must be altered by the back reaction 
of fundamental particles (including gravitational field) on the vacuum state; (iii) at very high tem- 
perature T ~ m p , thermal fluctuations are much more important than quantum-field fluctuations, 
and the vacuum is not stable. 



C. Back reactions on the vacuum 



The second part (OlAT^lO) in the total energy-momentum tensor 7^, (|8|) is the variation of the 
vacuum energy-momentum, 

(0|AT^|0} = (0|f^|0) s -(0|f^|0) ; (11) 

where (0|T^|0)s is the modified energy-momentum tensor caused by back reactions of external 
sources: (i) boundary conditions; (ii) classical fields; and (iii) fundamental particles on the vacuum. 

The vacuum energy-momentum is altered, due to external sources back interacting with vacuum 
fluctuations. The vacuum fluctuations, whose wavelengths are the order of the characteristic length 
scale of external sources, are most strongly affected. As a result, the variation of vacuum energy- 
momentum (0|AT^|0) is determined by the characteristic scales of external sources. In the case 
of external source being boundary condition, for example the Casimir effect [5|, the variation of 
vacuum energy-density is 

7T 2 1 

A ^Casimir _ 

Ae ° " ~720 

where the characteristic scale is Iq is the separation between two conducting plates. In the case 
of external source being classical gravitational field, for example the Hawking effect for 
an astrophysical black hole, the variation of vacuum energy-density is (a scalar field with two 
polarizations) 

Ae Hawking = -^Tfj, (13) 

where the Hawking temperature Th = l/(8irGM), M is the black hole mass and the characteristic 
scale is the black hole size GM. With macroscopic lengths lc > 10 _4 cm and astrophysical black hole 
sizes GM > 10 5 cm, one finds that the variations of vacuum energy-density are very much smaller 
than the energy-density m\ of an electron at rest, where m e is the electron mass. The variation of 
the vacuum energy-momentum (0|AT^|0) due to back reactions has very small gravitating effect 
and is negligible in the right-handed side of gravitational field equation ([7]). 
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In the assumption that vacuum fluctuations are homogeneously distributed in whole Universe, 
it is conceivable that the variations of vacuum energy-momentum induced by back reactions of 
external sources: topology boundaries of the Universe and classical gravitational perturbations at 
the cosmological scale (~ 10 28 cm) are extremely small, 

Ae Cosmo _ 10 -164[ Ge y]4_ ( 14 ) 

Comparing it with the energy-density e\ ([2j), one finds that Eq. (|14[) is of 117 order smaller. 
Therefore the variations of vacuum energy-momentum due to back reaction have entirely negligible 
impacts on the evolution of the Universe governed by gravitational field equation ([7]). 

In conclusion, the vacuum energy-momentum © (OlT^jO) = and the variation of the vacuum 
energy-momentum (jlip (0| AT^|0) ~ 0, the total energy-momentum tensor ([8]) becomes 

Tab T ab . (15) 

The gravitational field equation ([7]) can be written as the Einstein equation, 

Rab ~ ^9abR + ^9ab = 8irGT ab . (16) 

The cosmological term is clearly distinguished from the energy-momentum tensor of fundamental 
particles, and they play different roles in either side of the Einstein equation (|16p . It can be 
conceptually misleading, if one simply moves the A-term from the left-handed side to the right- 
handed side of the Einstein equation d7J) . 

D. Cosmological term ^ vacuum- energy term 

One is certainly allowed to optionally adopt a non-vanishing vacuum energy-momentum tensor 
(0|f^|0) + 0, and model the vacuum of quantum field theories as a perfect fluid with energy-density 
eo, pressure po and the equation of state eo = — Po, 

(0|^|0) = p g ab - (e +p )u a u b = -e g abl (17) 

where u a is the four velocity of fluid elements. However, in this framework, the energy-momentum 
tensor (0|T a {,|0} in Eq. (|l(jp is no longer the classical energy-momentum tensor T ab of fundamental 
particles, 



(0\f ab \0)=f ab ^T ab , 



(18) 
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and the total energy-momentum ©is 

Tab - (0| (f ab + |0) = T a6 - e <7a6, (19) 

where the back reaction term (0|AT^|0) is neglected. The Einstein equation becomes 

Rab ~ ^9abR + J^Qab = 8ttG (f ab - e g a b) ■ (20) 

It is completely not justified that one simply lets T ab = T ab , moves the cosmological term Ag ab to 
the right-handed side of the Einstein equation and relates it to the vacuum energy-density, 

Rab ~ \gabR = 8ttG (r ab - -^9ab^ , A = 8vrGe . (21) 

The problem of conceptually confusing the cosmological term and vacuum energy in the Einstein 
equation is much worst than the problem of large order magnitude discrepancy between the ob- 
served energy-density ([2]) and the Planck energy-density rut. 

In the empty space, there are no any fundamental particles T ab = 0, and the cosmological term 
is set to be zero, the vacuum Einstein equation (|16p reads 

Rab - \g ab R = o. (22) 

We attempt to study a possible origin of the cosmological term from quantum (thermal) fluctua- 
tions about gravitational instanton by using semi-classical saddle-point approach to the functional 
integral for quantum gravity. 



III. GRAVITATIONAL INSTANTONS AND FUNCTIONAL INTEGRALS 

The attractive nature of gravity that cannot be screened is the essential reason for many in- 
evitable instabilities of classical gravitating systems, such as the gravitational collapse and Jean 
instability. One might also worry about the instability of classical solutions of the Einstein equation 
against quantum fluctuations. To explore the structure of quantum (thermal) fluctuations about 

these classical solutions of the vacuum Einstein equation (I22D . one starts with the Euclidean action 

ft 

Se{9) and functional integral for quantum gravity [a], 

Z = Jv[g ab (x)}e- SEi9)+sauge - &xingteTms , (23) 
where the Euclidean action is given by, 

S E (g) = J (g) 1/2 drd 3 xR + boundary terms, (24) 
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where we make the Wick rotation, the Euclidean time dr = —idt and signature (+ + ++), 
(g^drd 3 x is the Euclidean volume element. Boundary terms are added to have an action which 
reproduces the vacuum Einstein equations under all variations of the metric that vanishes on the 
boundary. The functional integration is evaluated by integrating over all metrics g ab that are 
positive and definite and obey appropriate boundary conditions. The functional integral (|23H is 
treated by saddle-point methods. This is adequate for semi classical analysis of small quantum 
(thermal) fluctuations about classical solutions of the vacuum Einstein equation, in order to study 
the stability of these solutions against small quantum (thermal) fluctuations. 

The saddle-point evaluation of the functional integral Z (I23p starts by constructing the saddle- 
point of the action, namely, the classical solutions g s ab of the vacuum Einstein equation (|22p and 
R = 0. Expansion about these saddle points is performed by writing 

9ab{x) = g s ab (x) + 4>ab(x), (25) 

where <p a b are small quantum (thermal) fluctuation fields about g s ab as c-number background 
fields. The saddle-point metric g s ab (normally assumed to be nonsingular geodesically complete 
four-manifold) is colloquially termed a gravitational instanton [§]. Two types of boundary con- 
ditions that gravitational instantons g s ab (x) obey are discussed in Ref. [s]. One corresponds to 
the zero-temperature vacuum and another corresponds to the canonical ensemble at temperature 
T = 1/(3. 

The boundary conditions appropriate to the zero-temperature vacuum are termed by asymp- 
totically Euclidean (AE) [9J. An AE metric is one in which the metric approaches the flat metric 
on R 4 outside some compact set. Finite action requires the metric to be asymptotically, 

ds 2 = (l + ^ 5 ab dx a dx b + 0(r" 3 ), (26) 

where r is a four-dimensional radial coordinate and a a function of coordinates, but independent 
of r. The boundary at infinity is topologically S 3 . The gravitational instantons g s ab (x) obeyed such 
boundary condition (|26[) are colloquially termed gravitational AE instantons. 

The boundary conditions for the canonical ensemble at temperature T = 1/(5 are termed by 
asymptotically flat (AF) [9]. An AF metric is one in which the metric approaches the flat metric 
on R? x S 1 outside some compact set. Finite action requires the metric to be asymptotically, 



ds 2 = dr 2 + 1^1 + —j Sijdx'dx 1 +0( y \r\- i ); i,j = 1,2,3, (27) 

where |r| is a three-dimensional radial coordinate, a a function of coordinates, but independent of 
|r|; the Euclidean time r is a coordinate which is periodic with period (3. The boundary of infinity 
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is topologically S 2 x S 1 . The gravitational instantons g s ab {x) obeyed such boundary condition ([27]) 
are colloquially termed gravitational AF instantons. 

In Eqs. (|23|25p . treating <j) ab as small quantum (thermal) fluctuation fields and instantons g s ab 
as c-number background fields will generate the usual perturbation expansion. Up to the lowest 
order, i.e., the quadratic terms of quantum fields (p ab , one finds 

Z(g s ) « e ~ SE{9S) [vfiVfexp -~ j ' (g) 1/2 dTd 3 x^ M(g s )0 

= expj-W)-ilndet[-A% s )]}, (28) 

where metric indexes are omitted. SE(g s ) is the action of gravitational AE or AF instantons. The 
functional integration is evaluated by integrating over the quantum (thermal) fluctuation fields (p ab 
about the saddle-points g s ab . The matrix M{g s ) determines the properties of quantum fluctuation 
fields 4>ab about instanton g s ab . 

IV. QUANTUM (THERMAL) FLUCTUATIONS ABOUT FLAT SPACE 
A. The flat space at zero-temperature 

In the zero-temperature vacuum case, the positive-action theorem due to Schoen and Yau 
states that for any AE metric cf with R = 0, the action SE{g s ) is non-negative and SE(g s ) = if 
only if g^ b is flat. It was shown [8J that the action for any AE instanton must be zero, which follows 
from the fact that any AE instanton will be a solution of R ab = 0. The positive-action theorem then 
guarantees that AE instanton must be a flat metric. For zero-temperature, AE instanton is unique 
flat metric and one needs only examine the quantum fluctuations about flat space. The operator 
M{g s ) in Eq. (|28p is semi-positive definite 8j. Therefore the flat space in the zero-temperature 
vacuum case is stable quantum mechanically as well as classically. This precludes the possibility 
of flat space at zero temperature decaying by any mechanism. 

B. The flat space at non-zero temperature 

In the canonical ensemble at temperature T = 1/(3, the metrics g ab in Eq. (|23[) are strictly 
periodic in Euclidean time r with period (3: g a b(T, x) = g a b( T +f3, x), and integration over Euclidean 
time becomes, 

J dr = Np J dr, J (gf/Wx = N p j* dr J (g^tPx, (29) 
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and the Euclidean action ( 12 

Na 



S E (g) = -— J dr J (g) 1/2 d 3 xR + boundary terms, (30) 

where Na to be the number of periodic intervals in the Euclidean time. Defining the measure of 
functional integral as, 

N/3 

I V[g ab {x)} =UJ V[g ab {x)]p, (31) 
one rewrites the functional integral (|23)28p as the canonical partition function 

Z(g s ) = J Pb ab (x)]e- S ^)+gaugc-fixing 

ps ]Je~ SE ^ [[VtfVflpexp 



' ^dr ligYlWxtfMptf)* 



= Hexpt-S E (g s ) - ^lndet[^( 5 s )]| . (32) 

The classical solution g s , quantum (thermal) fluctuation field <j) and quadratic metric M.p{g s ) are 
strictly periodic in Euclidean time r with period (3. Therefore, the canonical partition function 
(|32|) is the same for different periodic sectors and one can normalizes Nr = 1. 

One defines the free energy J~{g s ) and the effective Euclidean action S E ^(g s ), to describe the 
equilibrium state of such a system, 

Z(g s ) = e^) = e -sgV) ; s f {gS) = ^ (/) + * hldeb [Mp(g')], (33) 

where S E ^(g s ) = —l3J-{g s ). We will see that the temperature T acts as an external " field " , and the 
real part and imaginary part of the effective Euclidean action describe the transition of space-time 
configurations. 

There is one periodic solution, namely, flat space g^ix) = 8 a b'. 

d 2 s = dr 2 + {dx 2 + dy 2 + dz 2 ). (34) 

This is a trivial AF instanton, since it has zero action SE{g s ) = 0. The semi-positive eigenvalues 
of the operator M(g s ) represent the modes of quantum (thermal) fluctuations about the classical 
solution g s ab {x) = 5 ab (134[) . which is a strictly local minimum of the action. The contribution of 
these modes to the free energy J-, to the lowest order in Eq. (|32p . will simply be the free energy 
of an ideal gas of relativistic massless particles with two helicities at temperature T, 
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the corresponding energy-density 



7T' 



P0 



rpl 

45 ' 



and the partition function 



cxp 



7T 

45" 



(36) 



(37) 



i.e., the contribution from thermal gravitons in a box of volume V at temperature T. The higher- 
order [in (167TG) 1 / 2 ] corrections {§] to the free energy, arising from the self-interactions for the gas 
of gravitons, leads to the infrared instability, i.e., Jeans instability. 



V. FINITE INSTANTON ACTION AT NON-ZERO TEMPERATURE 



Finite instanton action 



However, unlike the AE instanton in the zero-temperature vacuum case, the flat space fj34f) is 
not the unique AF instanton; there exist other periodic solutions of the Euclidean vacuum Einstein 



equations. A familiar instanton g^(x) is the Euclidean section of the Schwarzschid solution 11] : 



d 2 s = adr 2 + a l dr 2 + r 2 d8 2 + r 2 sin s 



a 



1 



2GM\ 1 / 2 
r J 



(38) 



in Boyer-Lindquist coordinates (t, r, 9, </>), where M is a mass parameter. The metric field must be 
periodic in r with period (3, as a result, M is not a fixed mass, but rather determined by inverse 
temperature (3 



P = 8irGM. 



(39) 



This gravitational instanton (|38j) is special classical field configuration characterized by its topo- 
logical properties: zero self-intersection number and an Euler number of two QjQ]- The classical 
Euclidean action of such an instanton is non-zero {§], 



Se(9 S ) = nPM, 



(40) 



at finite inverse temperature f3. This instanton action is not the contribution of a particle, or 
soliton, of mass M/2, since M depends on the temperature. 



12 1, which is the Euclidean 



The Schwarzschid solution (|38|) is a special case of the Kerr instanton 
section of the Kerr solution. It has been widely conjectured that the Kerr instanton is the unique 
AF instanton other than flat space (|34p []jj]. 
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The classical action of two instantons is given by 

4 2 V) = ^(/)- v z G X M y < 41 » 

where the first term is the action of two free instantons and the second term is their Coulomb 
interaction, and Az > 2GM is the distance between two instantons. This Coulomb interaction 
vanishes for a large separation of two instantons (Az 3> 2GM) and its maximum value for Az — > 
2GM is four times smaller than the action of two free instantons. 

Based on the assumptions that the spatial distribution of M instantons is homogeneous and the 
average distance between two instantons is {Az) > 2GM, up to the nearest neighborhood Coulomb 
interactions, the total classical Euclidean action is approximately given by, 

GM 



S^\g s )^AfS E (g" 



1 



(42) 



((Az) + 2GM)_ 

where MSE(g S ) is the classical action of M free instantons. In the following, we will neglect the 
Coulomb interaction and adopt the approximation, 

S { f\g s )*NS E (g s ), (43) 

since the correction due to the Coulomb interaction, i.e. the second term in Eq. (|42p . is less than 
25%. 



B. Vacuum to vacuum transition amplitude 



The functional integral Z (|23p represents the transition amplitude from vacuum to vacuum. The 
canonical partition function Z(g s ) (|32p is proportional to the transition amplitude from the vacuum 
of zero-instanton (5 a b hat space) to the vacuum of one-instanton g s at the finite temperature T. 
This transition amplitude is exponentially weighted by the finite classical action S E (g S ) of the 
instanton, 

z (l) {gS) = ffi,0|go6,0) T/ ^ = e -|/3M ) (44) 

[Oab, 0|O a fe,0) T=0 

where the amplitude is normalized by the amplitude in flat spacetime and at zero temperature [see 
Sect. (jIV Ap ] . Whereas the transition amplitude from the vacuum of zero-instanton to the vacuum 
of A/"-instanton sectors with different topological numbers is, 

(Sab,0\6 ab ,0) T =0 
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where the Coulomb interaction between instantons is neglected. The transition amplitudes (|44p and 
(|45|) indicates the configurations of instantons are weighted by e ~2@ M -^ in the functional integral 
([29]) and (l32l) . 

This is reminiscent of the t'Hooft gauge instanton. In the Euclidean space, the vacuum- 
to-vacuum amplitude of quantum gauge fields about the classical solution [15j is exponentially 



weighted by the action (8tt 2 /g 2 ) of the t'Hooft gauge instanton [l6], where g is a gauge coupling. 
This gives the transition amplitude of gauge instantons with different winding numbers. These 
studies have led to the 9- vacuum and resolution of the U(l) problem. 

When a thermal equilibrium is reached for a large Euclidean time interval At, and summing 
over all contributions of instanton sectors, the amplitude of the vacuum to vacuum transition (|45p 
has the form 



exp 



(A£ + i-)ArV 



(46) 



where AS (T) is a transition energy (rate) per unit of volume V. This relats to the real (imaginary) 
part of free energy T ([32]) or effective action Sj^(g s ) ([33]) . In next section, we will discuss that such 
vacuum transition is interpreted as the decay of the flat-space vacuum at a finite temperature, via 
the thermal nucleation of black holes. 



VI. INSTABILITY OF FLAT SPACE AT FINITE TEMPERATURE 

A. Instability of flat space 

The quantum (thermal) fluctuations 4>ab about the Schwarzschid instanton g s [see Eq. (1251) ] are 
determined by the eigenvalues and eigenfunctions of the operator [M.{g s )]j3 [see Eq. ([32]) ]. It was 



found 



13, 



18l ] that in addition to semi-positive eigenvalues, the operator A4(g ' )g has a negative 
eigenvalue A ~ — 0.19(GM) -2 , and corresponding mode (p^ b of quantum (thermal) fluctuation. 
The functional integral (|32p indicates that is an unstable mode for small quantum (thermal) 
fluctuations about the the Schwarzschid instanton g s . This instanton is therefore not a strict 
minimum, rather an unstable saddle point of the action. Its role in the thermodynamics of hot 
gravitons is similar to the top of the potential barrier. The rolling off the top of the potential 
barrier is caused by an infinitesimal perturbation. It was shown that the infinitesimal quantum 
(thermal) fluctuation ecf>^ b about the Schwarzschid instanton g , 

9ab{x) = gjj b (x) + afi b (x), (47) 
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decreases the Euclidean classical action Se(9 S ) (pTOJ) by [g] 

5S E (g S ) = -9.4 • W- 4 e 2 {GM)- 2 . (48) 

Consequently the functional integral (|32p is actually divergent, indicating the instability of running 
away from the classical action. 

This instability can be illustrated by Fig. (pQ) , where V is the energy potential of the Euclidean 
action. At finite temperature T, the Euclidean action has two extremes: (classical solutions of equa- 
tion of motion): one represented by A that is flat space solution 5 a b] another represented by B that 
the Schwarzschid instanton solution g~L. At a very low temperature T ^(^4) ~ 0, fluctuations 
about the state A of the flat space are not stable, due to both the classical Jean instability of long- 
wave length fluctuations and quantum mechanical tunneling of low-energy fluctuations through 
the potential barrier. While, at very high temperature, thermal fluctuations are dominate. The 
states at temperature T > V(B) degenerate to the state B of the Schwarzschid instanton g^L at 
temperature T = V(B), which simply sits at the top B of the potential barrier. These states are 
highly unstable, due to thermal fluctuations over the potential barrier decrease Euclidean action. 
As results, free energy decreases, and these states decays by gaining energy, i.e., rolling down from 
the top of energy potential. 

v 

B 



A 



FIG. 1: Instability of flat spacetime at finite temperature. V is the energy potential of Euclidean action. 
V{A) is the minimum of energy-potential and V(B) is the maximum of energy potential. 



B. Decay rate of flat space via thermal nucleation 

In Ref. [sj, this instability is interpreted as a decay of flat space (f34"|) at finite temperature T, 
and the decay process proceeds by quantum (thermal) fluctuations spontaneously nucleating black 
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holes of radius R = and mass M = (8-kGT)^ 1 . For all temperatures, the decay is only 

given by the saddle point at the top of the action, not given by the quantum tunneling process. 

This instability gives rise to an imaginary part in the free energy F{g s ) flU, which is related 
to the rate V (|46p of decay transition from the zero-instanton (flat) vacuum to the 7V-instanton 
vacuum. Since the decay process proceeds by nucleating black holes, this decay rate T is interpreted 
as the rate of thermal nucleating black holes per unit of volume, 



~ {g) l / 2 dTd?x' 

at a finite temperature T. In Refs. P, [if]], the thermal nucleation rate V (|49p is given by, 



(49) 



Im 



7T 



(50) 



where ojq is determined by the negative eigenvalue A, 

-A ~ 



1.74\ 2 



(51) 



for all temperatures. 

To obtain the thermal nucleation rate T, one has to calculate the functional integral Z f|23|) that 
integrates over the contributions from all instanton sectors J\f. In Ref. [§], the canonical partition 
function ([32]) for one instanton sector was calculated, 

M I 3 / 2 



212/45 



27T/3 



VZ u exp -S E (g b ) 



(52) 



where the imaginary factor of i/2 occurs from the one normalizable negative mode cj)\, Z° (|37p is 
the contributions from the stable quantum (thermal) fluctuations about the instanton, the factor 
arises from the renormalization counterterm and the regular mass /x should be taken 
of order the Planck mass m p . Neglecting the classical Coulomb interactions (|40p between the 
instantons, one estimated the contribution from the AMnstanton sector, 



Z^(g s ) 



V 



.64tt 3 G 3 / 2 

Summing over all instanton sectors, one obtained 

z(g s ) 

~ Z u exp 
and the free energy T = f3~ l In Z . 



M 



Z°exp \-NS E (sf 



(53) 



N=0 



2 WV U4vr 3 G 3 / 2 



cxp 



Se(9 S 



(54) 
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The thermal nucleation rate T per unit volume is thus given by Eqs. (|5QI54p . 



r-^(/,/3) 212 / 45 (^|)exp -S B tf) . (55) 



This formula gives the nucleation rate for pu re g ravity. It was extended to a general theory 
containing any number of massless matter field [20J , 



r.^a^-pf-J^V (56) 

v T ' 64vr 3 P ^ 16?rT 2 J v ' 

1 / 233 7 \ 

9 = — ( 212n 2 - —^3/2 - 13ni + -n 1/2 + n J , (57) 

where n s is the number of massless spin-s fields, the regular mass \i ~ m p and the classical action 
Se{9 S ) = ?Tip/(167rT 2 ). The thermal nucleation rate fj56|) reaches its maximum at high-temperature 
T ~ m p , and exponentially suppressed at low-temperature T <C m p . This indicates that thermal 
nucleation process could be very important in the very early Universe. 

Once these black holes have been nucleated, the unstable mode of small quantum (thermal) 
fluctuation e(j)\ will correspond to the subsequent expansion (or collapse) of the black hole as it 
absorbs (or emits) thermal radiation. Thus these black holes are colloquially termed quantum black 
holes in this article. The variations of the quantum black hole mass M and area A = IQir(GM) 2 

Q 

are given by [8|] 

M + SM, 5M = 9.4-10~ 3 e(GM)- 1 ; (58) 

A + SA, 5A = 9.4 • 10 -1 eG ! , (59) 

where positive and negative values of e ~ M 2 G correspond respectively expansion and collapse 
of the quantum black holes. These quantum black holes absorb (or emit) thermal radiation, then 
expand (or collapse), depending on the temperature of these quantum black holes T = l/(8irGM) 
being larger (or less) than the temperature of their environment. Henceforth, we adopt the notation 
QBH standing for quantum black hole. 

VII. GRAVITATIONAL INSTANTONS AND COSMOLOGICAL TERM 
A. The flat Robertson- Walker metric 

Suppose that the space-time symmetry of Universe is described by the Robertson- Walker line 
element with zero curvature k = and the scale factor a(r), 

d 2 s = dr 2 + a 2 (r)(dx 2 + dy 2 + dz 2 ) = a 2 (r?) \dr] 2 + dx 2 + dy 2 + dz 2 } , (60) 
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where the Euclidean time interval dr = —idt, and conformal time interval dtj = dr/air). The 
Robertson- Walker metric g^x) = a { T )$abi which is a space-time geometry that differs from the 
flat space-time geometry only by an overall scale factor a(r) > 0. 

Suppose that Universe is described by a canonical ensemble at temperature T = 1/(3, the 
Robertson- Walker metrics g^ b in Eq. (I60|) is strictly periodic in Euclidean time r with period (3: 
g^ b (r,x) = g^ b {r + /3,x). The metric g^ b (rj,x) = g^biv + ^-> x ), is strictly periodic in the Euclidean 
conformal time r\ with period 1, where we use the entropy conservation aT = 1 and temperature 
T is in unit of the Planck mass m p . 

The flat Robertson- Walker metrics g^ b (|60p is a classical solution to the vacuum Einstein equa- 
tion (}22|) and the Ricci scalar R = 0, provided, 

3 o/ da da 

- 2 a(r) = 0, a=- = a-. (61) 

The scale factor a is a time-independent constant, at a fixed temperature T. Thus the flat 
Robertson- Walker metrics g^ b = a5 a f), is a trivial AF instanton, which has zero action SE(g R ) = 0. 
The classical solution g R is a strictly local minimum of the action (^4 in Fig. {!])). The discussions 
on this local minimum can be found in section pVBp . 

B. Finite action of J\f instantons in functional integral 



As discussed in previous sections (jVUVip . there is another non-trivial AF Schwarzschid instan- 
ton (f38j) with the finite action SE(g S ) (0Q|) at temperature T. There is a topologically distinct 
configuration of AA-instanton sector whose finite action is approximately given by the finite action 
of M free instantons, i.e., J\fSE(g s ) [see Eq. (|42p]. The flat Robertson- Walker metrics g^ b = a5 a b at 
temperature T is not stable, against quantum (thermal) fluctuations. This flat Robertson- Walker 
spacetime is bound to spontaneously decay, via the thermal nucleation of QBHs. The rate of ther- 
mal nucleation process is given by Eq. (|56|) . This is an energy-gain process by decreasing finite 
action. 

Since the instanton action J\fSE(g S ) is finite and the thermal nucleation rate is nonzero, these 
topologically distinct configurations of instantons are not vanishing, and must be included into 
configurations of functional integral f|23|) and (|32|) [8|] . This is the same as the reason for including 



t'Hooft instantons into the functional integral in non-Abelian gauge theories 
integral (|32|) for quantum gravity at finite temperature T = /3 _1 becomes, 



2l|. Thus, functional 



Z({3) = f p[ 5(ifc ( x )] e -^5 E (sS)_ 5B ( 9 ) +gaU g e _ nxin gtermS ) ^ 

J instantons 
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where functional integral is carried over AMnstanton configurations. The Euclidean action Se(9) 
is given in Eq. (I30p . This indicates that the effective action for quantum gravitational field at 
finite temperature contains contributions from gravitational instantons. 

As discussed in section ()Vip , that the finite action of Af-instanton sector in functional integral 
(1621) indicates the decay process by thermal nucleating QBHs. The decay rate T at a finite tem- 
perature T is given by Eqs. (|49|55|56p . independently from space x and time r. Therefore, the 
total number of QBHs nucleated M is given by 

TV = TV At, (63) 

where At is the Euclidean time interval when thermal nucleation process takes place. The finite 
action of AMnstantons in the canonical partition function (|62p is then, 

MS E {g S ) = TV At, (64) 

where the finite action for a single instanton Se(9 S ) = 0M/2 (|38p is time-independent. 

C. Partition function for time-evolution of canonical ensemble 

So far we have discussed the functional integral (partition function) describing a canonical 
ensemble at a fixed temperature. Now we turn to study the functional integral describing the 
adiabatic evolution of a canonical ensemble, where temperature T = is a function of time. 
We consider a sequence of canonical ensembles S(tq), S(t±), Sfa), S(ts), ■ ■ -,S(tn), at different 
temperatures T(tq) > T(t{) > Tfo) > T{t%) > ■ ■ ■ > T(rjv), and scaling factors gs(to) < a(ri) < 
0(^2) < o(r3) < • • ■ < a(rjv), at different times To < t\ < T2 < T3 < • • • < tn, in comoving frame. 
In this case, both the mass parameter M(t{) (|39p and thermal nucleation rate T(rj) (|56p of QBHs 
are functions of time T{. At each instantaneous moment r, and temperature T{ri), the canonical 
partition function is given by the functional integral (|29I3QI32[) and (|62p with normalization Ng = 1. 
For an infinitesimal time interval dr = — t%, we assume that the time-evolution of canonical 
ensemble S(t{) is strictly adiabatic and total entropy of canonical ensemble is conserved, i.e., 
a(Ti) = (3(Ti). 

We define the functional integral for the time-evolution of canonical ensembles as a time-ordering 
product of canonical partition functions (|62|) for canonical ensembles at different temperature 

T( Ti ) = p-Hn), 

N 

z = rn*)] = z[p( To )] ■ z[(3( n )] ■ zipfo)] ■ ■ ■ z[p( TN )}. (65) 

i=0 
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Further, we define the integration over Euclidean time sequence t% as, 

'TO i=Q Jb i = Q 

four-dimensional volume element as, 

i=0 i=0 ' ' 

and the measure of functional integral over different time period (3 fa) as, 

N P{n) 



rr N rP{n) N 

dr'^Y. dr' = Y;Pfa), r^r N , (66) 

J to j=0 JO i=Q 

element as, 

J (<?) V V* = £ flr*)^) = J2J Q ' dr'j ' G?) 1 / 2 ^*, (67) 



i=0r'=0 

Then, the functional integral (|65p can be written as, 

N P{n) 



J V[g ab (x)} = f[ f[ / V[g ab (x,T>)}. (68) 



z= n n f^[gabfa,T')]e- MsE ^ s) - sEi9)+gauee ^ ine , 

i=0 t'=0 



(69) 



where the Euclidean action Se(9) and finite instanton action MSE(g S ) are respectively given by 
the sum over their counterparts at different temperature Tfa). The Euclidean action Se^) is, 

1 N rP{n) r 

S E (g) = -Jq^q E / dT ' J (9) 1/2 d 3 xR + boundary terms. (70) 



i=0 



The finite TV- instanton action MSE(g S ) ([Ml) becomes, 

(71) 



r 1 



A Ti -Mfa)Tfa 



NS E {g s )=Y d Pfai)Vfa i 

i=0 

The factor J2iLo Pfa)Vfa) gives the four dimensional volume (f6"7|) in the function integral (|68I69|) . 
Because the mass parameter M and decay rate V in Eq. ([7T]) are time-dependent, we should 
integrate Mfa)Tfa)/2 over the time-interval At, = Tj — To when thermal nucleating process occurs, 

Ar±Mfa)Tfa) - jT dr^M^Or^'). (72) 

As results, the finite AMnstanton action J\fSE{g s ) (fTT]) is given by 

A^(/) = | {g) l l 2 d A xp h fa,r), (73) 

where 

PA(r ,r)= [ dT' l -M(r')T{T ! ). (74) 

J TO ^ 

Actually, Pa(t ,t) is the energy-density of AA QBHs nucleated in the Euclidean time interval 



19 



D. Minkowski formulation 

Performing the back Wick rotation from the Euclidean signature (+ + ++) to Minkowski one 
( — h ++), we can obtain the effective action and functional integral in the Minkowski formulation. 
The Euclidean action (|70p becomes the Minkowski one, 

S(g) = ^-77 / {-g) l/2 d A xR + boundary terms, (75) 

167rG J 

where {—gY^d^x is the Minkowski volume element, and functional integral measure (j68|) is changed 
correspondingly. The finite action of A^-instanton actions (|T3|) becomes the cosmological term, we 
write it as 

Sa(9 S ) = J(-g) l/2 d 4 x PA (t ,t), (76) 



where 



and r(t') is the thermal nucleation rate per unit of Minkowski volume [see Eq. (j39 

r _ M 

{- g y/ 2 dtd?x' 

The cosmological constant is then given by, 



PA(to,t) = f dt'\M{t')T{t'), (77) 

Jto 1 



(78) 



A = 8irGp A (t ,t). (79) 

The above formulations and discussions can be straightforwardly extended to the case that 
matter fields are present. The total effective action in functional integral is then, 

S(g) + S A (g s ) + S M , (80) 

where Sm is the action for fundamental particles, whose energy momentum tensor is T a ^. In fact, 
the non-zero temperature T is attributed to the thermal energy of relativistic particles. Based on 
the principle of least action, we obtain the Einstein equation ([7]) with the cosmological term given 
by Eq. (|76p . Such cosmological term has two features: (i) its geometric origin from the nucleation 
of gravitational instantons, and appearance in the right-handed side of Einstein equation ([7]); (ii) 
its energetic origin from the energy-gain of nucleating QBHs 5E A = —5E, where 5E < is the 
vacuum-energy variation of the flat space at finite temperature. This temperature is due to thermal 
energy of relativistic particles. The energy-momentum tensor corresponding the cosmological term 
is T^ b = —g a bPA- The total energy conservation is given by (T a b + T^ b )' ,a = in the process of 
nucleating QBHs. 



20 



VIII. EARLY EVOLUTION OF UNIVERSE 
A. Inflation 



We consider the post-Planckian era of very early Universe, at the time t > to = 1, the scale 
factor a(io) > a (to) = 1 an d the temperature T < Tq = 1, where and henceforth, the temperature 
T and time t are in the Planck unit. The temperature is due to the thermal graviton gas with 
energy-density po (i36j) . Suppose that the early Universe is given by the spatial-flat geometry of 
the Robertson- Walker line element with the scale factor a{t) and zero curvature k = (I60h . and 
the Einstein equation ([7|) describing Universe expansion becomes, 

H2 ^(S) =^-iPo(t)+PA(t)]. (81) 
where ph{t) is given by Eq. (|77[) with the initial time to = 1 and temperature Tq = 1. 



PA® = \ f dTM{T)T[T{r% 



0.87mt rt ( \ 



W2^J lo dr {W)) e (82) 
which is the energy-density of QBHs nucleated from initial time to an d temperature To to final 
time t and temperature T(t). Corresponding to the QBHs' energy-density (f82"j) . the QBHs' energy- 
density is given by, 

rt 

The Universe expands and its temperature T(t) decreases, a(t)T(t) = 1 for the entropy- 
conservation. As results, the energy-density po(t) (|36|) decreases, whereas pA(t) dBZD increases 
and asymptotically approaches a constant p\ for i > 1. Thus, Eq. (j81j) implies an inflationary 
Universe. How does such inflation phase end? 



n A (t) = \ f dTY[T[r)\, 

2 J to 



B. Hawking radiation and reheating. 

When the Universe temperature T(t) is smaller than the temperature T(r) = 1/ {8itGM{t)) of 
QBHs that are created at an earlier time r < t, i.e., T(t) < T(r), these QBHs lose their masses 
by the Hawking radiation. On the other hand, accretion occurs if T{t) > T(r) (see Eqs. (|58l) 
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and (|59p ). If we were only to consider emission and absorption of gravitons, the mass- variation of 



QBHs is given by 



6, 
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^ = ^[T\t)-T\T)]A.R\Tl (84) 

where the black hole size R(t) = 2M(r). For T(r) 3> T(t), we approximately obtain, 

M H (t) ~ M(r)[l - ^T 3 (r)(t - r)} l '\ (85) 

indicating that QBHs lose their mass-energy. We speculate that the extremal minimum of QBHs' 
mass-energy is the order of the Planck scale. 

The Hawking process (I85p occurs, contemporaneously with QBH nucleation. Due to the Hawk- 
ing radiation, the energy-density p A (t) ([82]) of QBHs is reduced. The variation of the energy-density 
p\(t) can be obtained by Eq. ([82]) . 

ApA(i) = \ f drAMT[T{ T )], (86) 
2Jt 

0.87mt ft / 1 \ i 

= f / drAM — - e «^ 2 M, 

128vr 3 J t0 vT(r)y 

where 

AM = M# (t) - M(r), (87) 

and Mn(t) is given by Eq. (|85p . The energy-density variation of QBHs Ap/^(t) < 0. The energy- 
density pa(£) in Eq. (f8Tj) should be replaced by 

p(t)=p A (t) + Ap A (t) = \ f dTM H (r)T[T(r)}. (88) 

This shows that the energy-gain p\(t) in QBH nucleation is converted into the radiation energy 
\ApA(t)\, which starts to reheat the Universe. 

Assuming a(t)T(t) = 1 and 9 = 212/45 (gravitons), to = 1 and To = 1, we use the method 
of iterating procedure in the time "t" to numerically integrate Eqs. (|81l88p for the evolution of 
Universe, by taking into account both the nucleation process (j56!82j) and the Hawking process 
(|85l86p . At the beginning (t < O(10 2 )), the nucleation process dominates over the Hawking 
process. The variation of energy-density p(t) (I88p is rather slow, so that the "graceful exit" problem 
is avoided and the evolution of Universe is inflationary. However, at the end (t > O(10 2 )), the 
variation of energy-density p(t) (|88p is rather rapid, the Hawking process is dominant, eventually 
ends the inflation and reheats Universe. 
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FIG. 2: The energy-densities p m (t), pA(t), p(t) and |A/9a(^)I as functions of time. 

C. Numerical results. 

In the initial phase t ~ 1, the energy-density of QBHs p\{t) ([82]) is negligible, compared with 
the energy-density of thermal graviton gas po (136p , The solution to Eq. (|8ip is radiative, 

ai=a t 1/2 , T^Tot- 1 / 2 , (89) 

where ao ~ 1 and To ~ 1 are initial scaling-factor and temperature. a^To ~ 1 implies that the 
initial entropy So = (aoTb) 3 is given by 0(1) quantum states of Planck energy in the Planck volume. 
As the time t increases, po(t) decreases, p\(t) increases and becomes dominant in Eq. (f8T1) . 

In Fig. [2], we plot po(t) (cross line) and p\(t) (short-dash line) as functions of the time t. It 
is shown that (i) the pre-inflationary phase (|89p for t < 10; (ii) the inflationary phase for t > 10, 
where po (j36|) is vanishing and p\(t) is approaching to an asymptotic value 

p A s 9.03 • 10" 2 m^, (90) 

and the corresponding QBHs number-density (|83p is ~ 10 _2 mp. In Fig. [2], we also plot the 
radiation energy |Apa(£)| ([86]) (dot-line) and the energy-density p(t) ([88]) (long-dash line), showing 
that p{t) slowly varies for 8 < t < 100 and p(t) — > for t > 100. Correspondingly, in Fig. [3], we 
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plot a(t) and T(t) ~ l/a(t), which show that (i) the pre-inflationary phase (|89p for t < 4, (ii) an 
exponential inflation a(t) ~ aoexp(iV e £) for 4 < i < 110 and (iii) a{i) approaching 10 30 for t > 110. 

We consider that the inflation ends at t ~ tf, when the QBHs' mass-energy (I85|88p is reduced 
to Mu — l/(8vr) corresponding to the QBH temperature Tu ~ 1. We find that tf ~ 113, af ~ 
2.84 • 10 28 and the e-folding factor N e t f ~ 65.5. 

The reheating process occurs when t ~ tf. We assume that apart from the Hawking radiation 
from each QBH, these QBHs collide each other and annihilate into other particles, since the energy- 
density p{t) and number-density n\ of QBHs is still very high, comparable with Planck densities, at 
the end of inflation phase t ~ tf. We define that the reheating phase ends at time th, when t — > th, 
QBHs' densities p(t) — > and n\(t) — > 0. In the other wards, the reheating phase ends at t = th, 
when the QBHs' energy-density p\ (|90p has been completely converted into the energy-density of 
other relativistic particles, and QBHs has completely annihilated into other relativistic particles. 
The reheating temperature Th at the end of reheating t = th can be estimated by the thermal 
energy-density of thermal relativistic particles, 

pZ = g s ^T* = p A ~\Ap A (t h )\, (91) 

where g s stands for the summation over contributions of all relativistic particles created in the 
Hawking process and QBHs annihilation. Since all possible relativistic particles are created, the 
reheating temperature must be smaller than the initial temperature, i.e., Th < Tq ~ 1. We are 
not able to determine Th and th, and leave them as parameters in this article. For Th ~ O(10~ 2 ), 
an enormous entropy Sh = (T^a^) 3 ~ O(10 79 ) is produced. In the time interval (tf,th) for the 
reheating phase, \Ap A (t)\ — > p\, the scaling factor a(t) slowly varies, see Fig. [3]. 

These numerical results Figs. |2|3| depend on the initial temperature To and time to. In partic- 
ular, the initial temperature To plays a crucial role, for its exponential dependence exp — l/(16-7rT 2 ) 
in QBH nucleating rate (|82p . While, numerical results are not sensitive to the parameter 9 ([3]). 



IX. COSMOLOGICAL CONSTANT AT THE PRESENT UNIVERSE. 

After reheating to the temperature Th, the nucleation of QBHs starts again, and the density of 
energy-gain p\(t) is given by Eq. (f82j) with the initial time th and temperature Th- 



p\{t) = \f drM(r)T[T(r)], 
2 Jt h 



0.87mt ft f 1 



1024vr 4 J th \T(t 



dr —— e ie - T2 M, (92) 
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FIG. 3: The scaling factor a(t) and temperature T(t) — l/a(t) as functions of time. 

where temperature T is smaller than the reheating temperature Th (T < Th). Because the reheating 
temperature Th is lower than the initial temperature To (Th < To ~ 1), QBHs' density p\(t) (f92j) 
is exponentially suppressed by lower temperature T < 1, thus smaller than p$ (|91|) the thermal 
energy-density of thermal relativistic particles, 

PA«Po^7t t>t h . (93) 

Therefore, the thermal energy-density pfi dominates in Eq. (|8ip for the evolution of the Universe. 
As a consequence, the Universe begins the evolution described by the Standard Cosmology: 

a(t) = a h (t/t h ) a , T(t)=T h (t/t h y a , (94) 

with total entropy Sh = (a^T^) 3 . The index a is determined by radiation or matter dominate 
epochs. We consider this as a new era initiated with th,Th and ah, independently from the infla- 
tionary era before the reheating. 

The energy-density p\(t) is mainly contributed from QBHs nucleated in the reheating. Due to 
the Hawking radiation ([85 1). the variation of energy-density Ap^(t) is given by Eq. (j86p with initial 
time th and temperature Th, 

A/4(i) = \f drAMT[T(r)\, (95) 



25 



and Ap A (t) < 0. Analogously to Eq. ([88]) . the total energy-density of QBHs nucleated after the 
reheating ear is, 

Phit) = p h A (t) + Ap h A (t) 

0.87m* ft ( 1 N^" 1 ) i , 

f / dTM H (——) e le-r-W . (96) 



128.3 A, " VT(r) 
This is related to the "dark-energy" density, and the cosmological constant is given by 

A = 8nGp h (t). (97) 

Assuming that the mass of QBHs has been reduced to the minimal mass Mh = 1/(8.) at the 
present time t 3> th, we obtain, 

/N 0.87mt, ft ( 1 V - 1 ) i , 

from Eq. (fMj) . Substituting solutions a(i) and T(t) (f94j) into /%(£) (|98j) . we obtain, 

/N 0.87m* r* ( T * \ ^ 

Ph (t) = £ / dr e ""Vh 

m; 1024. 4 4 ^T^y 

^(1^ 1/a 
= 1024^ 2a fe r ( ff ^^) (") 

where the incomplete Gamma-function is given by 

T(8,z 2 ,z x ) = f ^ dxx s - l e~ x , (100) 

and 

S = i^tl, (101) 
t 2a 1 

z 2 = TTT^^Z , Zl 



The asymptotic representation of the incomplete Gamma- function (jlOOp , 

r(a,32,«i) ^^ _1 e- zi . (102) 

for Z2 S> -Zi S> 1 and ^2 - > oo. Using this asymptotic representation, we obtain the "dark-energy" 
density and cosmological constant at the present time t S> th, 

0.87 mfeh (3-0) ~Ti^5f A Q ~ nn „s 
^ - 64 ^3 2a T h ' e S A = 8vrGp ft . (103) 

which are approximately independently of time t. Analogously to Eq. (|83f) . we compute the number 
density of QBHs created after the reheating phase, 

0.87ml r t / i 

drr[T(T)j = 



n h (t) = / drr[T(r)] = / t dr [f^ ) e T5TO ' (J0I) 
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and obtain the number-density of QBHs at present time, 



3+, ._ L 

16tvT: 



0.87 mph ( 3 - 



We find that the "dark-energy" density Ph(t) (f9"8|) . cosmological constant (j97|) and the number- 
density of QBHs increase, when the Universe is in radiation domination (t > t^,T < T/J, and 
asymptotically approaches the constants (I103P and fjlQ5|) . when the Universe is in matter dom- 
ination (t S> ift,T <C r^). The numerical values of these constants (|103j) and (|105j) depend on 
three parameters Th,th and 0, in particular, crucially and sensitively depends on the reheating 
temperature Th via the exponential function exp — l/(167rT^). 

With the reheating time th — 113, the present time t ~ 10 61 , the index a = 1/2 in the radiation 
dominant phase, and 6 = 203/45 for the particle content of the Standard Model, we obtain 

p h ~ 6.7 • 10" 120 m^; A ~ 1.68 • lO -118 ™ 2 ,, (106) 

by setting reheating temperature Th — 8.45 • 10 -3 . This is self-consistent within this theoretical 
model and agrees with present observations (|1|2|) . The corresponding residual QBHs number- 
density is 

n h ~ 2.67 • 10" 22 /cm 3 , (107) 

which is rather small so that the annihilation probability of residual QBHs can be neglected at the 
present time. The total number N of residual QBHs is an approximate constant in the present 
evolution of Universe. 



X. EQUATION OF STATE 

As Universe adiabatically expands, the internal energy E of particles decreases, following the 
energy-conservation, 

P 5V + 5E = 0, (108) 

where p is the pressure, 5V > is the variation of Universe volume and 5E < the variation of 
internal energy of particles. The Equation of State is 

V = (7 - IK (109) 

where e is the energy-density in the co-moving frame and thermal index 7 > 1. Being different 
from particles moving upon the manifold of four-dimension space-time, residual QBHs are geomet- 
rical holes embed in the manifold. When the manifold is stretched in the expansion of Universe, 
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these geometric holes are stretched and their sizes become larger, as we will see, their energies 
increase, rather than decrease. This implies that the "Equation of State" of residual QBHs must 
be completely different from the equation of state (11Q9H of particles. 

The radius and mass of residual QBHs are R = 2M and M = l/(8irT) (T ~ m p ). Their 
3-dimensional surface and volume are A = AttR 2 and V = 47ri? 3 /3. As the manifold is stretched in 
the expansion of Universe, the radius R, the surface A and volume V of QBHs are stretched to be 
larger. The volume variation of a residual QBH is, 

SV = A5R = 2A5M > 0, (110) 

where 5R and SM are the variations of residual QBHs' radius and mass. The thermal energy 
(kinetic energy) of residual QBHs can be neglected, compared with their masses at the Planck 
scale. Eq. (IllOp shows that 5R > and 5M > for 5V > 0, indicating residual QBHs increase 
their size and mass-energy in the expansion of Universe. The energy-conservation for a residual 
QBH is, 

P 5V + 5M = 0, (111) 

in the adiabatic expansion of Universe. In order to obtain " Equation of State " of residual QBHs, 
we assume that 

• residual QBHs are uniformly distributed in the entire manifold of Universe; 

• the manifold is uniformly stretched in the expansion of Universe. 

Therefore, in the expansion of universe, the energy variation of Af residual QBHs is 8{MM) = M5M 
and the volume variation of residual QBHs is 8(NV) = AfSV. The energy-density of residual QBHs 
is then given by, 

6(NM) 5M 1 . , 

where the thermal energy (kinetic energy) of residual QBHs' is neglected. From Eqs. (jllllll2p . 
we obtain the " Equation of State " of residual QBHs, 

p = -e. (113) 

Taking into account the small correction due to the thermal energy (kinetic energy) of residual 
QBHs, we can write the Equation of State (|113p as p = — (1 — c)e, where c is the radio of residual 
QBHs' thermal energy and mass-energy, and c < 1. The small parameter c must be positive, 
since the thermal energy (kinetic energy) e& of residual QBHs decreases, i.e., 5ek < 0, as Universe 
volume expands 5V > 0. 
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XI. SUMMARY AND SOME REMARKS 

In an empty space and at zero temperature, classical particles' energy-momentum tensor (jlOD 
T a b = 0, and has no contribution to the vacuum Einstein equation (122j) . We have argued in Sec. 
(|II Bp that quantum fields' vacuum energy-momentum ([9]) (0|TX|0) = 0, and has not contribution 
to the source term (r.h.s.) of Einstein equation (|16p . This is in accordance with the positive-energy 
theorem [10], which actually requires zero mass-energy in the Minkowski spacetime, and the ground 
state is the flat Minkowski spacetime. This implies that the vacuum energy of quantum fields is 
not gravitating and thus the cosmological term is very unlikely related to the vacuum energy of 
quantum fields. 

In this article, we study the possibility that the cosmological term is given by the energy-density 
of QBHs, which is related to the finite action of classical gravitational instantons. The nucleation 
of QBHs is attributed to the unstable quantum fluctuations about the classical gravitational in- 
stantons, i.e., the decay of flat space at finite temperature. Based on such cosmological term, we 
study the evolution of the early Universe and present value of the cosmological constant. These 
are still very preliminary results, and there are many questions for further studies in future. 
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